Method and apparatus for determining receiver beam forming vectors in multi-user MIMO systems

ABSTRACT

According to one embodiment there is provided a method of determining receiver beam forming vectors for a plurality of receivers in a MIMO system. Each receiver comprises a plurality of receiver antennae. The MIMO system also comprises a plurality of transmitter antennae. The method comprises quantifying the properties of the channels of the MIMO system in a channel property matrix including interference between individual channels and performing Cholesky decomposition of the receiver property matrix row-by-row into a lower triangular matrix with unit elements along its diagonal, a diagonal matrix and the Hermitian transpose of said triangular matrix. Performing the Cholesky decomposition comprises determining, when performing a step of the Cholesky decomposition for a row, a receiver-beam forming vector associated with the row that maximises the component of the diagonal matrix in the row.

CROSS-REFERENCE TO RELATED APPLICATION(S)

This application is based upon and claims the benefit of priority from UK Patent Application No. 0911901.7, filed on 8 Jul. 2009, the entire contents of which is incorporated herein by reference.

FIELD

Embodiments described herein relate generally to a method and apparatus for use in a multi-user multiple input multiple output (MIMO) system. Embodiments described herein relate to a method and an apparatus for determining one or more of receiver beam forming and transmitter precoding and feedback matrices for use in such a system.

BACKGROUND

Multi-user multiple input multiple output precoding techniques are attractive for practical MIMO system deployment thanks to their potential to increase the system's capacity. In multiple input multiple output systems employing precoding, channel knowledge is used at the transmitter in order to enhance the link quality. A recent development relates to nonlinear precoding/multi-user MIMO. This technique shows promise for improving the system capacity. Among nonlinear precoding schemes, Tomlinson-Harashima precoding (THP) is one of popular schemes due to its simplicity and good performance. A potentially practical scheme is the minimum means square error (MMSE) THP proposed by K. Kusume, M. Joham, W. Utschick, and G. Bauch in “Cholesky factorization with symmetric permutation applied to detecting and precoding spatially multiplexed data streams,” IEEE Trans. on Signal Proc., vol. 55, no. 6, pp. 3089-3103, 2007, the entirety of which is incorporated herein by this reference. This scheme provides good performance in terms of bit error rate (BER).

BRIEF DESCRIPTION OF THE DRAWINGS

An embodiment of the present invention will now be described by way of example only and with reference to the accompanying drawings, in which:

FIG. 1 shows a block diagram of a multi-user Tomlinson-Harashima precoding (THP) system;

FIG. 2 shows a sequence diagram of an embodiment;

FIG. 3 shows a pseudo code for one part of the Cholesky decomposition method of an embodiment;

FIG. 4 shows pseudo code for a user selection algorithm;

FIGS. 5 and 6 show simulation results comparing a system with receiver beam forming vectors determined based on the embodiment with known systems; and

FIG. 7 shows simulation results comparing a system with receiver beam forming vectors determined based on and users selected using the embodiment with known systems.

DETAILED DESCRIPTION

According to one embodiment there is provided a method of determining receiver beam forming vectors for a plurality of receivers in a MIMO system. Each receiver comprises a plurality of receiver antennae. The MIMO system also comprises a plurality of transmitter antennae. The method comprises quantifying the properties of the channels of the MIMO system in a channel property matrix including interference between individual channels and performing Cholesky decomposition of the receiver property matrix row-by-row into a lower triangular matrix with unit elements along its diagonal, a diagonal matrix and the Hermitian transpose of said triangular matrix. Performing the Cholesky decomposition comprises determining, when performing a step of the Cholesky decomposition for a row, a receiver-beam forming vector associated with the row that maximises the component of the diagonal matrix in the row.

The channel property matrix may also comprise an indication or quantification of the noise expected to be present in the channels. The plurality of receivers may be non-cooperating receivers. The receiver beam forming vectors of all of the receivers of the MIMO system may thus be expressed in a block diagonal matrix.

The component of the diagonal matrix may be maximised by determining the left singular vector of a matrix quantifying the characteristics of the channel associated with a said element of the diagonal matrix. The left singular vector may be the left singular vector associated with the maximum singular value of the matrix. The matrix quantifying the characteristics of the said channel may quantify channel gain, channel noise and the said channel's interference with channels of receivers for which a beam forming vector had previously been determined.

According to another embodiment there is provided a method of determining a receiver beam forming vector for a receiver of a plurality of receivers in a MIMO system, wherein each receiver comprises plural receiver antennae. The method comprises forming a matrix quantifying the properties of a channel of the receiver and the interference the said receiver experiences from other receivers for which a receiver beam forming vector had previously been determined. The receiver beam forming vector can then be determined. This includes determining the left singular vector of the matrix that quantifies the properties of a channel of the receiver and the interference the said receiver experiences from other receivers.

The matrix quantifying the characteristics of the said channel may be:

${H_{i}H_{i}^{H}} + {\gamma\; I_{N_{r}}} - {\sum\limits_{j = 1}^{i - 1}{s_{i\; j}s_{i\; j}^{H}q_{j}}}$ for the second and subsequent rows/steps of the Cholesky decomposition, wherein H_(i) is the channel matrix of the receiver associated with the current row/step of the Cholesky decomposition, H_(i) ^(H) is the Hermitian transpose of H_(i), γ is the expectation of the noise that will be received, I_(N) _(r) is an N_(r)×N_(r) identity matrix, wherein N_(r) is the number of receiver antennae of the receiver, s_(ij) is

$s_{i\; j} = \left\{ \begin{matrix} \frac{H_{i}{\overset{\Cap}{h}}_{j}^{H}}{q_{1}} & {if} & {j = 1} \\ \frac{\left( {{H_{i}{\overset{\Cap}{h}}_{j}^{H}} - {\sum\limits_{n = 1}^{j - 1}{s_{in}l_{j\; n}^{*}q_{n}}}} \right)}{q_{j}} & {if} & {j > 1} \end{matrix} \right.$ ĥ_(j) ^(H) is the Hermitian transpose of r_(j)H_(j), wherein r_(j) is the receiver beam forming vector associated with the j-th element of the diagonal matrix, H_(j) is the channel matrix associated with the j-th element of the diagonal matrix, q_(J) is:

$q_{i} = {\Psi_{i\; i} - {\sum\limits_{k}^{i - 1}{l_{i\; k}l_{i\; k}^{*}q_{k}}}}$ wherein:

${l_{i\; j} = {\frac{1}{q_{j}}\left( {\Psi_{i\; j} - {\sum\limits_{n}^{j - 1}{l_{in}l_{j\; n}^{*}q_{n}}}} \right)}},{{{for}\mspace{14mu} i} > j}$ l_(i j) = 1, for  i = j l_(i j) = 0, for  i < j $\Psi_{i\; i} = {{{{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{i}^{H}} + {\gamma\mspace{14mu}{and}\mspace{14mu}\Psi_{i\; j}}} = {{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{j}^{H}}}$ ${{\overset{\Cap}{h}}_{i} = {r_{i}H_{i}}},$ s_(ij) ^(H) is the Hermitian transpose of s_(ij), l_(jn)* is the conjugate of l_(jn). The j-th element may be an element that has previously been determined. N_(r) may be the number of receiver antennae of each receiver. For the second and subsequent rows/steps of the Cholesky decomposition the sum term in the matrix quantifying the characteristics of the channel may be zero as interference with other channels needs to be taken into account. The receiver beam forming vector may have unit norm.

The matrices provided by the Cholesky decomposition of the channel property matrix Ψ can be expressed as:

$\begin{matrix} {\Psi = {LQL}^{H}} \\ {= {\begin{pmatrix} 1 & 0 & \cdots & 0 \\ l_{21} & 1 & \ldots & 0 \\ \cdots & \cdots & \cdots & \cdots \\ l_{K\; 1} & l_{K\; 2} & \cdots & 1 \end{pmatrix}\begin{pmatrix} q_{1} & 0 & \cdots & 0 \\ 0 & q_{3} & \cdots & 0 \\ \cdots & \cdots & \cdots & \cdots \\ 0 & 0 & \cdots & q_{K} \end{pmatrix}\begin{pmatrix} 1 & l_{21}^{*} & \cdots & l_{K\; 1}^{*} \\ 0 & 1 & \cdots & l_{K\; 2}^{*} \\ \cdots & \cdots & \cdots & \cdots \\ 0 & 0 & \cdots & 1 \end{pmatrix}}} \end{matrix}$ wherein L is the lower triangular matrix with unit elements along its diagonal, Q is the diagonal matrix and L^(H) is the Hermitian transpose of L, wherein q₁ to q_(k) are the elements of the diagonal matrix defined as:

$q_{i} = {\Psi_{i\; i} - {\sum\limits_{k}^{i - 1}{l_{i\; k}l_{i\; k}^{*}q_{k}}}}$ wherein:

${l_{i\; j} = {\frac{1}{q_{j}}\left( {\Psi_{i\; j} - {\sum\limits_{n}^{j - 1}{l_{in}l_{j\; n}^{*}q_{n}}}} \right)}},{{{for}\mspace{14mu} i} > j}$ $\Psi_{i\; i} = {{{{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{i}^{H}} + {\gamma\mspace{14mu}{and}\mspace{14mu}\Psi_{i\; j}}} = {{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{j}^{H}}}$ ${\overset{\Cap}{h}}_{i} = {r_{i}H_{i}}$ wherein l_(ik)* is the conjugate of l_(ik), H_(i) is the channel matrix of the receiver associated with the current row/step of the Cholesky decomposition, γ is the expectation of the noise that will be received and r_(i) is the receiver beam forming vector associated with the current row/step of the Cholesky decomposition. I_(ij) may be determined after each determination of q_(i).

The above method is such that it allows determination of the receiver beam forming vectors for all desired receivers before precoding and/or feedback matrices that may be used in the transmitter have to be determined. Such a transmitter may be a transmitter that employs Tomlinson-Harashima precoding (THP), for example Minimum Means Square Error (MMSE) THP. The method may thus further comprise determining a precoding and/or a feedback matrix for a transmitter using Tomlinson-Harashima precoding after two or more of the receiver beam forming vectors have been determined. The precoding and/or feedback matrices may be determined after all of the receiver beam forming vectors have been calculated.

The method may further comprise receiving signals at one or more of the plurality of antennae of a said receiver and combining the signals using the receiver beam forming vector.

The above discussion has focussed on scenarios where the number of users participating in the MIMO system may correspond to the number of users that can participate in the system. Frequently, however, the number of users attempting to participate in a MIMO system exceeds the maximum number of users of such a system.

In such as situation it may be necessary to exclude users from participation in data transmission in the MIMO system. It is desirable that only users maximising transmission quality are selected for data transmission in the MIMO system. As discussed above, the selection of a correct receiver beam forming vector can minimise the mean square error associated with signal transmission. The receiver beam forming vector that minimises MSE is the receiver beam forming vector that maximises an associated component of the diagonal matrix calculated by Cholesky decomposition. It would thus be desirable to determine from among a number of thus far unselected receivers the one receiver that maximises the associated element of the diagonal matrix created by Cholesky decomposition of a channel property matrix quantifying not only the channel properties associated with a receiver in question but also the properties of a channels associated with previously selected receivers.

To achieve the above, a channel property matrix may be defined for two or more or even for each thus far unselected receiver. This channel property matrix may quantify the properties of the channel associated the unselected receiver as well as the properties of the channels of all thus far selected receivers. Cholesky decomposition may then be performed for the channel property matrix and a receiver beam forming vector maximising the element of the diagonal matrix associated with the unselected receiver can be determined. This determination also provides a value for the associated element of the diagonal matrix. As this computation can be performed for two or more or even for each of the thus far unselected receivers, the receiver providing the maximum element of the diagonal matrix can be determined. This receiver is the one of the thus far unselected receivers under consideration that minimises the mean square transmission error. This receiver may then be selected.

According to another embodiment there is provided a method of selecting a user from a plurality of un-selected users for participation in signal transmission in a MIMO system comprising a plurality of transmitter antennae and a plurality of receivers. Each receiver in the system comprises a plurality of receiver antennae. The method comprises performing individually for more than one of the unselected users a quantification of the properties of the receiver channel of the unselected user and of receiver channels of previously selected users, if any, in a receiver property matrix. This quantification includes a quantification of interference between individual channels. Cholesky decomposition of the receiver property matrix row-by-row into a lower triangular matrix with unit elements along its diagonal, a diagonal matrix and the Hermitian transpose of said triangular matrix is then performed for each of the individually determined receiver property matrices. When a step of the Cholesky decomposition is performed for a row, a receiver beam forming vector associated with the row is determined such that the receiver beam forming vector maximises the component of the diagonal matrix in the row. The method further comprising selecting from among the more than one unselected users the user having the maximum diagonal matrix component. Selection of users for participation in data transmission in the MIMO system provides a further opportunity to improve performance.

Each of the thus far unselected users may be considered in this manner, although doing so is not essential and only a smaller sub-set of the thus far unselected users may instead be considered.

Any of the above methods may further comprise receiving signals at the plurality of antennae of a said receiver and combining the signals using the receiver beam forming vector.

According to another embodiment there is provided a device for determining receiver beam forming vectors for a plurality of receivers in a MIMO system in which each receiver comprises a plurality of receiver antennae and that further comprises a plurality of transmitter antennae. The device comprises a quantifier for quantifying in a channel property matrix the properties of receiver channels of the MIMO system, including interference between individual channels and a processor arranged to perform Cholesky decomposition of the received matrix in a row-by-row fashion into a lower triangular matrix with unit elements along its diagonal, a diagonal matrix and the Hermitian transpose of the lower triangular matrix. The processor is further arranged to determine, when performing a step of the Cholesky decomposition for a row, a receiver beam forming vector associated with the row that maximises the component of the diagonal matrix in the row.

Embodiments also extend to a transmitter for use in a MIMO system, wherein the transmitter comprises apparatus for selecting a user or users from a plurality of un-selected users. Embodiments also extend to a MIMO system comprising such a transmitter.

According to another embodiment there is provided a device for determining a receiver beam forming vector for a receiver of a plurality of receivers in a MIMO system, wherein each receiver comprises plural receiver antennae. The device comprises a quantifier for quantifying in a matrix the characteristics of a channel of the receiver and the interference the said receiver experiences from other receivers for which a receiver beam forming vector had previously been determined. The device further comprises a processor arranged to determine the receiver beam forming vector, wherein the processor is arranged to determine the left singular vector of the matrix.

Embodiments also extends to a transmitter for use in a MIMO system, wherein the transmitter comprises the device for determining the receiver beam forming vectors. Embodiments also extends to a MIMO system comprising the device for determining the receiver beam forming vectors.

According to another embodiment there is provided an apparatus for selecting a user from a plurality of un-selected users for participation in signal transmission in a MIMO system comprising a plurality of transmitter antennae and a plurality of receivers, wherein each receivers comprises a plurality of receiver antennae. The apparatus is arranged to, individually for more than one of the unselected users, quantify the properties of the receiver channel of the unselected user and of receiver channels of previously selected users, if any, in a receiver property matrix including interference between individual channels. The apparatus is further arranged to, again individually for more than one of the unselected users, perform Cholesky decomposition of the receiver property matrix row-by-row into a lower triangular matrix with unit elements along its diagonal, a diagonal matrix and the Hermitian transpose of said triangular matrix. When a step of the Cholesky decomposition is performed for a row a receiver beam forming vector that is associated with the row is determined, wherein the receiver beam forming vector is determined such that it maximises the component of the diagonal matrix in the row. The apparatus is further arranged to select from among the more than one unselected users the user having the maximum diagonal matrix component.

FIG. 1 shows a block diagram of a multi-user Tomlinson-Harashima precoding (THP) system. The system comprises a transmitter 10 with K transmitter antennae in communication with K receivers 20. Each receiver 20 has N_(r)≧1 receiver antennae. The channel matrix characterising the channels between all of the transmitter antennae and the receivers 20 can be expressed as: H=[H ₁ ^(T) , . . . , H _(K) ^(T)]^(T) εC ^(KN) ^(T) ^(×K)  (1) wherein H_(k)εC^(N) ^(T) ^(×K)  (2) is the k^(th) receiver's channel. These receiver channels of the individual receivers are indicated in FIG. 1 by reference numerals 30.

The transmitter 10 of the THP system generates a spatial stream vεC^(K×1) from a data stream vector dεC^(K×1) that is to be transmitted to the receivers 20. K is precoded using precoding matrix FεC^(K×K) (identified by reference numeral 40) to form signals x, for transmission to the receivers 20 via the channels 30. The signal y received at the receivers 20 can be modelled as: y=HFv+n  (3) wherein n is the noise vector for all the receivers: n=[n ₁ ^(T) , . . . , n _(K) ^(T)]^(T) εC ^(KN) ^(T) ^(×1)  (4)

The elements of n have zero-mean Gaussian distribution.

The spatial streams v is generated from the signal vector d through use of a feedback filter 60 according to: v=B⁻¹d  (5) where BεC^(K×K) is a lower triangular matrix with unit diagonal elements.

The modulo operation 60 shown in FIG. 1 is introduced to reduce the transmitter power increase due to the implementation of the inverse of B. Equation (5) thus becomes: v=B ⁻¹(d+τa)  (6) wherein a is a complex integer vector and τ is a positive real number depending on the modulation schemes. For QAM modulation τ could be: τ=2(|c| _(max)+Δ/2)  (7) where |c|_(max) is the absolute value of the real or imaginary part of the modulation constellation symbol with greatest magnitude, and Δ is the smallest distance between two modulation constellation symbols. The vector of a can be found by a simple modulo operation and the detail of the modulo operation can be found in section IV of the above incorporated paper by Kusume.

Turning now to the receiver, the data stream at each receiver could be estimated by equalization as: {circumflex over (d)}=β ⁻¹ Ry=β ⁻¹ R(HFv+n)  (8) where β is a scalar used to normalize the transmit power. In THP the receivers are assumed to act independently from each other. For this reason the receiver beam forming matrix RεC^(K×K·N) ^(r) is a unity matrix for N_(r)=1, i.e. if each receiver only comprises a single receiver antenna. If each receiver comprises more than one receiver antennae, that is for N_(r)>1, the beam forming matrix R is a block diagonal matrix:

$\begin{matrix} {R = \begin{bmatrix} r_{1} & 0 & 0 & 0 \\ 0 & r_{2} & 0 & 0 \\ 0 & 0 & \cdots & 0 \\ 0 & 0 & 0 & r_{K} \end{bmatrix}} & (9) \end{matrix}$ with r_(k)ε

^(1×N) ^(T) is the receiver beam forming vector for k^(th) receiver.

The error between the received signal {circumflex over (d)} and the data stream vector d that was transmitted and the covariance matrix of this error are respectively: ε={circumflex over (d)}−d  (10) Φ_(εε)=E└εε^(H)┘  (11) where E[.] denotes the expectation.

The total mean square error (MSE) can be obtained as: ρ=E(tr(εε^(H)))=tr(E(εε^(H)))=tr(Φ_(εε))  (12) where tr (•) represents trace of matrix. The optimal R, F, B and β with minimum mean square error (MMSE) criteria can be obtained by minimizing the mean square error ρ and the problem formulation to find the optimal R, F, B and β is as follows:

$\begin{matrix} {\left\{ {F_{opt},R_{opt},B_{opt},\beta_{opt}} \right\} = {\underset{F,R,B,\beta}{argmin}\left\{ {{tr}\left( \Phi_{ɛɛ} \right)} \right\}}} & (13) \end{matrix}$ and subject to the following condition: tr(FΦ _(vv) F ^(H))=E _(S)  (14) where: Φ_(vv)=E└vv^(H)┘  (15) and E_(S) is the transmit power.

Equation (8) can be rewritten as: {circumflex over (d)}=β ⁻¹ Ry=β ⁻¹ R(HFv+n)=β⁻¹(ĤFv+{circumflex over (n)})  (16) wherein: Ĥ=RH=[(r ₁ H ₁)^(T), . . . , (r _(K) H _(K))^(T)]^(T)=[(ĥ ₁)^(T), . . . , (ĥ _(K))^(T)]^(T) and  (17) {circumflex over (n)}=Rn

The noise power is kept unchanged by constraining the receiver beam forming vector so that it has unit norm. This keep the covariance matrix of {circumflex over (n)} that same as that of n, i.e. Φ_({dot over (n)}{dot over (n)})=Φ_(nn).

Kusume et. al., referred to above, provides a MMSE solution for the case that each receiver comprises a single receiver antenna only, that is for N_(r)=1. As mentioned above, in this case the receiver beam forming matrix is simply the identity matrix and does not need to be determined. Kusume characterises the transmission of signals by defining: Ψ=(ĤĤ ^(H) +γI)=LQL ^(H)  (18) where L is a lower triangular matrix with unit diagonal and Q is a diagonal matrix. L and Q can be jointly calculated using a Cholesky decomposition. Moreover:

$\begin{matrix} {\gamma = \frac{{tr}\left( \Phi_{\overset{\Cap}{n}\overset{\Cap}{n}} \right)}{E_{S}}} & (19) \end{matrix}$ that is the expectation of the noise received together with the signal after combination of the received signals using the receiver combination vectors and normalised by the transmitter power E_(S). In a practical system, the covariance matrix Φ_({dot over (n)}{dot over (n)}) of {circumflex over (n)} can be determined once a transmitter has acquired knowledge of the noise power of the receivers. γ can then be computed using equation (19).

The feedback matrix and the precoding matrix can be obtained according to:

$\begin{matrix} {\overset{\_}{B} = L} & (20) \\ {{\overset{\_}{F} = {{\overset{\_}{\beta}\;\hat{F}} = {\overset{\_}{\beta}\;{\overset{\Cap}{H}}^{H}L^{- H}Q^{- 1}}}}{wherein}} & (21) \\ {{{\hat{F} = {{\overset{\Cap}{H}}^{H}L^{- H}Q^{- 1}}}\overset{\_}{\beta} = \sqrt{\frac{E_{S}}{{{\hat{F}\left( {\text{:},1} \right)}}_{2}^{2} + {\sigma_{v}^{2}{{\hat{F}\left( {\text{:},{2:K}} \right)}}_{F}^{2}}}}}{{{where}\mspace{14mu}\sigma_{v}^{2}} = \frac{\mu}{\mu - 1}}} & (22) \end{matrix}$ represents the precoding loss due to modulo operation at receiver and μ is the modulation order.

In a MIMO system comprising receivers with one receive antenna each the receiver beam forming matrix is the identity matrix, as discussed above. Based on F, B and β determined using equations (20) to (22), the MSE of equation (12) can thus be calculated as:

$\begin{matrix} {{\rho\left( {\overset{\_}{R},\overset{\_}{F},\overset{\_}{B},\overset{\_}{\beta}} \right)} = {\gamma{\sum\limits_{i = 1}^{K}\frac{\sigma_{v\; j}^{2}}{q_{i}}}}} & (23) \end{matrix}$ where σ_(vj) ²=1 for j=1 and σ_(vj) ²=σ_(v) ² for j>1 and q_(i) is the i-th diagonal element of matrix Q in equation (18).

The above pre-coding matrix and feedback matrix are generated with given receiver ordering. Kusume et. al., referred to above, also discloses two receiver ordering algorithms for further improving the performance of the disclosed method. These algorithms are summarised in tables VI and VII of Kusume et. al. respectively, which are incorporated herein by this reference.

Techniques for generating the pre-coding and feedback matrices exist. These techniques, however, make assumptions that relax the constraints in equation (13).

As mentioned above, the MMSE solution suggested by Kusume et. al. focuses on a scenario where each receiver comprises a single antenna, N_(r)=1. In this case the receiver beam forming matrix R is the identity matrix and thus does not need to be separately determined. If the receivers 20, however, comprise more than one receiver antennae, N_(r)>1, then it is necessary for the receiver beam forming matrix R to be determined, before F, B and β can be obtained based on equations (20), (21) and (22). As R, F, B and β are mutually embedded in equations (20), (21) and (22), F, B and β cannot be obtained using these equations.

METHOD OF DETERMINING RECEIVER BEAM FORMING VECTORS

In the following a method is described in which the receiver beam forming matrix R is determined based on Cholesky decomposition. In each step of this method the receiver beam forming vector r₁ . . . r_(K) for a given receiver is generated by minimising the mean square error in equation (22).

Equation (17) can be re-written as:

$\begin{matrix} \begin{matrix} {\Psi = \left( {{\overset{\Cap}{H}{\overset{\Cap}{H}}^{H}} + {\gamma\; I}} \right)} \\ {= {LQL}^{H}} \\ {= {\begin{pmatrix} 1 & 0 & \cdots & 0 \\ l_{21} & 1 & \ldots & 0 \\ \cdots & \cdots & \cdots & \cdots \\ l_{K\; 1} & l_{K\; 2} & \cdots & 1 \end{pmatrix}\begin{pmatrix} q_{1} & 0 & \cdots & 0 \\ 0 & q_{3} & \cdots & 0 \\ \cdots & \cdots & \cdots & \cdots \\ 0 & 0 & \cdots & q_{K} \end{pmatrix}\begin{pmatrix} 1 & l_{21}^{*} & \cdots & l_{K\; 1}^{*} \\ 0 & 1 & \cdots & l_{K\; 2}^{*} \\ \cdots & \cdots & \cdots & \cdots \\ 0 & 0 & \cdots & 1 \end{pmatrix}}} \end{matrix} & (24) \end{matrix}$

Based on Cholesky decomposition q_(i) and l_(ij) can be obtained as:

$\begin{matrix} {q_{i} = {\Psi_{i\; i} - {\sum\limits_{k}^{i - 1}{l_{i\; k}l_{i\; k}^{*}q_{k}}}}} & (25) \\ {{l_{i\; j} = {\frac{1}{q_{j}}\left( {\Psi_{i\; j} - {\sum\limits_{n}^{j - 1}{l_{in}l_{j\; n}^{*}q_{n}}}} \right)}},{{{for}\mspace{14mu} i} > j}} & (26) \end{matrix}$ wherein: Ψ_(ii)=ĥ _(i) ĥ _(i) ^(H)+γ  (27) Ψ_(ij)=ĥ_(i)ĥ_(j) ^(H)  (28)

As can be seen from equations (25) and (26), q_(i) is successively generated step by step and is affected by the parameters of previous receivers and of the current receiver. During each step to generate q_(i) the MSE of equation (23) is minimised by selecting r_(i) so as to maximise q_(i).

It will be appreciated from the above that successive determination of the receiver beam forming vectors is possible as Cholesky decomposition of the channel property matrix provides a lower triangular matrix with unit elements along its diagonal. Receivers associated with a particular row of the channel property matrix consequently experience interference from signals to receivers in subsequent rows. However, each receiver interferes with receivers associated with subsequent rows of the channel property matrix. The receiver beam forming vectors of a receiver can therefore be determined taking only the properties of the receivers in preceding rows into account. A row-by-row determination of the receiver beam forming vectors is thus possible. The method may be used for determining receiver beam forming vectors for receivers in a data transmission system using Tomlinson-Harashima precoding, MMSE-based linear precoding or MMSE-based vector perturbation precoding, to name a few.

It will be appreciated that, as a consequence of the above, the receiver associated with the first row of the channel property matrix, or with the first element of the diagonal matrix respectively, can be considered in isolation from the other receivers. The channel of this receiver, however, interferes with other channels and these interferences are taken into account when determining subsequent receiver beam forming vectors. By maximising the elements of the diagonal matrix the means square error of the received signals is minimised.

As a consequence of the manner in which the channel property matrix is decomposed, the only variable in performing the decomposition of a row of the matrix is the receiver-beam forming vector. The receiver-beam forming vector can therefore be chosen to minimise the transmission error. This can be achieved by maximising the element of the diagonal matrix that is being determined by the current step of the Cholesky decomposition routine (as mentioned above). The elements of the associated row of the lower triangular matrix that is also determined during the step of the Cholesky decomposition process can be determined after the element of the diagonal matrix has been determined. A receiver beam forming vector for a receiver can thus be determined based solely with knowledge of the properties of the associated receiver channel, of receiver channels associated with receivers for which receiver beam forming vectors had previously been determined and of interference between these channels. The use of a zero forcing constraint is not required in the above method. The method of the embodiment may increase spatial diversity, giving rise to better system performance.

FIG. 2 shows a flow chart illustrating how this can be achieved. In a first step 100 the number K of receivers 20 is defined and γ is determined. The channel matrix H is defined according to equations (1) and (2) above.

It can be seen from equation (24) that for the first receiver q₁=Ψ₁₁ applies. Ψ is defined above in equation (24) and according to this definition Ψ_(ii) =ĥ ₁ ĥ ₁ ^(H)+γ  (29) applies. Equation (29) can further be re-written as: Ψ_(ii) =ĥ ₁ ĥ ₁ ^(H) +γ=r ₁(H ₁ H ₁ ^(H) +γI _(N) _(r) )r ₁ ^(H) =r ₁Γ₁ r ₁ ^(H)  (30)

All of the components of Γ₁ have been determined in step 100 shown in FIG. 2. The Hermitian transpose of the receiver beam forming vector r₁ of the first receiver that maximises Ψ₁₁ and therefore q₁ can thus be determined in step 110 of FIG. 2 as the dominant left singular vector of Γ₁. Once the receiver beam forming vector r₁ of the first receiver has been determined, ĥ₁ is calculated based on: ĥ₁=r₁H₁  (31)

As can be seen from equation (24) and in analogy with equation (30) above, q₂ is:

$\begin{matrix} \begin{matrix} {q_{2} = {\Psi_{22} - {l_{21}l_{21}^{*}q_{1}}}} \\ {= {{{\overset{\Cap}{h}}_{2}{\overset{\Cap}{h}}_{2}^{H}} + \gamma - \frac{{\overset{\Cap}{h}}_{2}{\overset{\Cap}{h}}_{1}^{H}{\overset{\Cap}{h}}_{1}{\overset{\Cap}{h}}_{2}^{H}}{q_{1}}}} \\ {= {{r_{2}\left( {{H_{2}H_{2}^{H}} + {\gamma\; I_{N_{r}}} - \frac{H_{2}{\overset{\Cap}{h}}_{1}^{H}{\overset{\Cap}{h}}_{1}H_{2}^{H}}{q_{1}}} \right)}r_{2}^{H}}} \\ {= {r_{2}\Gamma_{2}r_{2}^{H}}} \end{matrix} & (32) \end{matrix}$

The components of Γ₂ are again known from the outputs of the preceding steps, steps 100 and 110. In step 120 the Hermitian transpose of the receiver beam forming vector r₂ for the second receiver can therefore be determined by calculating the dominant left singular vector of Γ₂, that is the left singular vector associated with the maximal singular value of Γ₂. Based on r₂ it is possible to calculate ĥ₂ in analogy with equation (17) and l₂₁ based on equation (26).

As can be seen from equation (24) and in analogy with equations (30) and (32) above, q₃ can be calculated as:

$\begin{matrix} \begin{matrix} {\quad{q_{3} = {\Psi_{33} - {l_{31}l_{31}^{*}q_{1}} - {l_{32}l_{32}^{*}q_{2}}}}} \\ {= {{{\overset{\Cap}{h}}_{3}{\overset{\Cap}{h}}_{3}^{H}} + \gamma - \frac{\Psi_{31}\Psi_{31}^{*}}{q_{1}} - \frac{\left( {\Psi_{32} - {l_{31}l_{21}^{*}q_{1}}} \right)\left( {\Psi_{32} - {l_{31}l_{21}^{*}q_{1}}} \right)^{H}}{q_{2}}}} \\ {= {{{\overset{\Cap}{h}}_{3}{\overset{\Cap}{h}}_{3}^{H}} + \gamma - \frac{{\overset{\Cap}{h}}_{3}{\overset{\Cap}{h}}_{1}^{H}{\overset{\Cap}{h}}_{1}{\overset{\Cap}{h}}_{3}^{H}}{q_{1}} -}} \\ {\frac{\left( {{\overset{\Cap}{h}}_{3}{\overset{\Cap}{h}}_{2}^{H}{\overset{\Cap}{- h}}_{3}{\overset{\Cap}{h}}_{1}^{H}l_{21}^{*}} \right)\left( {{\overset{\Cap}{h}}_{3}{\overset{\Cap}{h}}_{2}^{H}{\overset{\Cap}{- h}}_{3}{\overset{\Cap}{h}}_{1}^{H}l_{21}^{*}} \right)^{H}}{q_{2}}} \\ {= {r_{3}\left( {{H_{3}H_{3}^{H}} + {\gamma\; I_{N_{r}}} - \frac{H_{3}{\overset{\Cap}{h}}_{1}^{H}{\overset{\Cap}{h}}_{1}H_{3}^{H}}{q_{1}} -} \right.}} \\ {\left. {\ldots\frac{\left( {{H_{3}{\overset{\Cap}{h}}_{2}^{H}} - {H_{3}{\overset{\Cap}{h}}_{1}^{H}l_{21}^{*}}} \right)\left( {{H_{3}{\overset{\Cap}{h}}_{2}^{H}} - {H_{3}{\overset{\Cap}{h}}_{1}^{H}l_{21}^{*}}} \right)^{H}}{q_{2}}} \right)r_{3}^{H}} \\ {= {r_{3}\Gamma_{3}r_{3}^{H}}} \end{matrix} & (33) \end{matrix}$

The components of Γ₃ are again known from the values output by the preceding steps, steps 100, 110 and 120. The Hermitian transpose of the receiver beam forming vector r₃ for the third receiver can therefore be determined by calculating the dominant left singular vector of Γ₃. Based on r₃ it is possible to calculate ĥ₃ in analogy with equation (17) and l₃₁ and l₃₂ based on equation (26).

From the above description and in light of equations (25) and (26) above, the person skilled in the art will appreciate that each receiver beam forming vector r_(i) can be calculated based on the values output by computation steps performed to calculate preceding receiver beam forming vectors r₁ to r_(i−1). In particular, the matrix Γ_(k) describing the properties of signal transmission to the receiver k in question is determined in the manner described for the first three receivers above. As will be appreciated from the above, this matrix also accounts for the interaction of the received signal with signals transmitted to the other receivers. The Hermitian transpose of the receiver beam forming vector r_(k) for receiver k is then determined by computing the dominant left singular vector of the matrix Γ_(k).

Based on the newly determined receiver beam forming vector r_(k) the matrix ĥ_(k) describing the properties of the channel and receiver beam forming matrix as all as l_(ij) per equation (26) can then be calculated, so that the receiver beam forming vectors r_(k+1) of the next receiver can then be determined.

FIG. 2 illustrates this by the loop operations indicated by reference numerals 130 to 150. In first step the number of the receiver k for which the receiver beam forming vector r_(k) is being determined is incremented in step 130. In step 140 a determination is made if k is smaller or equal the total number of receivers K as defined in step 100, that is if there are still some receiver beam forming vectors to be determined. If this is the case the procedure moves on to step 150, where the matrix Γ_(k) describing the properties of signal transmission to the receiver k in question is determined in the above described manner and based on equations (25) and (26) and the output values of the preceding computational steps. The Hermitian transpose of the receiver beam forming vector r_(k) for receiver k is then also determined in step 150 by computing the dominant left singular vector of the matrix Γ_(k). Based on r_(k) the matrix ĥ_(k) for the receiver in question as well as the terms l_(k,j) for j<k are then determined. The procedure then moves on to step 130 to initiate the computation of the receiver beam forming vector for the next receiver. If it is determined in step 140 that all receiver beam forming vectors have already been determined, then the determined vectors are output in step 160.

The determination of the values l_(ij) is based on equation (26). The computation of l_(ij) can be performed prior to determining a receiver beam forming vector r_(i) in a loop according to the following pseudo code shown in FIG. 3 for the second and higher number receivers, that is for r_(i) with i≧2. The loop is initiated in Step 200 and operates a number of times that corresponds to the number of receivers for which receiver beam forming vectors had previously been determined. The number of times the instructions in the loop are repeated thus corresponds to the number of entries in the relevant row of matrix L as shown in equation (24) above.

It will be appreciated that for any of the values l_(i,l), that is for the values in the first column of the matrix L, the sum term shown in equation (26) will be zero. Computation of the first value in each row of L from the second row onwards can therefore be simplified when compared to the computation of values in columns of L other than the first column. For this reason as a first step, step 210, within the loop it is determined whether or not the next value computed is in the first column of the matrix L. If this is the case, then computation is based on a simplified version of equation (26), as set out in step 220 shown in FIG. 3. All other values in the matrix L, that is all values in columns of this matrix other than the first column are computed using the entirety of equation (26), as indicated in Step 230 shown in FIG. 3.

Once all values l_(ij) have been determined, these values can be inserted into equation (25) and the matrix Γ_(i) required for the determination of the beam forming vector of the receiver that is currently being considered can then be derived according to:

$\begin{matrix} \begin{matrix} {q_{i} = {{{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{i}^{H}} + \gamma - {\sum\limits_{j = 1}^{i - 1}{l_{ij}l_{ij}^{*}q_{j}}}}} \\ {= {{r_{i}H_{i}H_{i}^{H}r_{i}^{H}} + \gamma - {\sum\limits_{j = 1}^{i - 1}{r_{i}s_{ij}s_{ij}^{H}r_{i}^{H}q_{j}}}}} \\ {= {{r_{i}\left( {{H_{i}H_{i}^{H}} + {\gamma\; I_{N_{r}}} - {\sum\limits_{j = 1}^{i - 1}{s_{ij}s_{ij}^{H}q_{j}}}} \right)}r_{i}^{H}}} \\ {= {r_{i}\Gamma_{i}r_{i}^{H}}} \end{matrix} & (34) \end{matrix}$

As discussed above, it is then possible to compute the Hermitian transpose of the receiver beam forming vector r_(k) ^(H) by calculating the dominant left singular vector of Γ_(i). Based on this output ĥ_(i) and l_(ij), i>j can be calculated in turn.

It will be appreciated from the above that q_(i) at the i^(th) step is always obtained as: q_(i)=r_(i)Γ_(i)r_(i) ^(H)  (35) where Γ_(i) is the function of parameters obtained from previous receivers and the channel of current receiver: Γ_(i) =f(H _(i) ,H _(j) ,q _(j) ,l _(k,n) ,r _(j)) j,k<i, n<j  (36)

The minimum mean square error in equation (23) is minimised if q_(i) is the dominant singular value of Γ_(i) and r_(i) ^(H) is the corresponding dominant singular vector.

METHOD OF SELECTING A USER/ORDERING USERS

The above described method determines receiver beam forming vectors wherein the receivers are placed in a predetermined order. It may, however, be desirable to alter the order in which the receiver beam forming vectors are determined. It may also be desirable for the most suitable user to be selected from among a plurality of users, for example in a situation where the number of user N attempting to transmit via the MIMO system exceeds the number of receivers K that is for N≧K. A receiver can hereby be called a user without loss of generality. In this case not all users can participate in the network and the most suitable users for participation in the network are thus preferably selected. For the purpose of selection the K users are given indices {π(1), . . . , π(K)} and the N users are given indices {1, 2, . . . , N}.

During the step i of the Cholesky decomposition process, a user π(i) is then selected from among the set of users (π(1), . . . , π(i−1))^(c) that remain unselected at the start of this step i. The user selected in this arrangement is the user that has the maximum q(i) from amongst the unselected users.

The pseudo code for this selection process is shown in FIG. 4. C^(c) is the complement set of C and C∪C^(c)={1, 2, . . . , N}. In step 300 the set C is emptied and a loop is entered in Step 310. This loop is repeated K times for the selection of K users from among the N potential users. In Step 320 a subset k comprising those of the N potential users that have thus far not been selected is defined. It is then determined in Step 330 whether a user has already been selected. If this is not the case, that is if i=1, then the values q_(i) ^(k) can be determined using equation (30). A value q_(i) ^(k) is determined for each of the remaining unselected users in subset k in Step 340. If the second user is being selected the values q_(i) ^(k) can be determined using equation (32) in Step 350. Again a value q_(i) ^(k) is determined for each of the remaining unselected users in subset k in Step 350. The selection of any user beyond the second user can be based on equation (34) in Step 370. In this step a value q_(i) ^(k) is determined for each of the remaining unselected users in subset k.

Once a value q_(i) ^(k) has been determined for each of the remaining users in the subset k (be that through Step 340, 360 or 370), the remaining unselected user with the maximum value q_(i) ^(k) is selected as the i^(th) selected user in Step 380. In step 390 the selected user π(i) is then added to the set C. As C^(c) is/remains the complement set of set C, this step of course reduces the number of remaining users available for selection in the next repetition of the loop shown in FIG. 4. Once this loop has been repeated K-times a full set of users C={π(1), . . . , π(K)} has been selected.

It will be appreciate that, because the users are selected using equations (32) and (33), the receiver beam forming vectors for each selected user needs to be calculated before further users can be selected. Thus, once the set of users have been selected all receiver beam forming vectors (with the possible exception of the receiver beam forming vector for the last selected receiver/user) have already been calculated.

Based on the selected set of users C and on the associated receiver beam forming vectors the Cholesky decomposition of ψ can be performed in accordance with equation (18) above, thereby providing the feedback matrix B and the pre-coding matrix F as well as β, using equations (20) to (22) respectively.

As discussed above, a row in a matrix that is decomposed/factorised using Cholesky decomposition does not interfere with the Cholesky decomposition of preceding rows. The receiver property matrices may thus be arranged so that the channel properties of an unselected user currently under investigation are placed in the matrice's last row. This has the advantage that Cholesky decomposition for the rows other than the last row does not have to be repeated. If previously selected users were selected by performing Cholesky decomposition of their receiver property matrices, then the result of this computation are already known and do not have to be determined again. These results can thus be re-used when determining the element of the diagonal matrix associated with the unselected user currently under investigation. It will be appreciated that this can lead to a considerable saving in expenditure of computational resources, in particular if the number of unselected users under investigation is high. The order of the users in the channel property matrix may thus be the order in which the users have been selected.

In an alternative arrangement the precoding matrix F and the feedback matrix B can be determined based on a re-ordered set of users as described in the following. Further details of this ordering method can also be found in Table VI of the Kusume et. al. paper incorporated by reference in its entirety above.

Let PΨ′P ^(T) =P(ĤĤ ^(H) +γI)⁻¹ =P ^(T) =L′ ^(H) Q′L′  (37) where P is a permutation matrix and

$\begin{matrix} {\overset{\Cap}{H} = {RH}} & (38) \\ {R = \begin{bmatrix} r_{\pi{(1)}} & 0 & 0 & 0 \\ 0 & r_{\pi{(2)}} & 0 & 0 \\ 0 & 0 & \ldots & 0 \\ 0 & 0 & 0 & r_{\pi{(K)}} \end{bmatrix}} & (39) \\ {H = \left\lbrack {H_{\pi{(1)}}^{T},H_{\pi{(2)}}^{T},\ldots\mspace{14mu},H_{\pi{(K)}}^{T}} \right\rbrack^{T}} & (40) \end{matrix}$ Then the precoding matrix F and feed-back matrix B can be obtained as:

$\begin{matrix} {B^{\prime} = L^{\prime - 1}} & (41) \\ {F^{\prime} = {{\beta^{\prime}\hat{F}} = {\beta^{\prime}{\overset{\Cap}{H}}^{H}P^{T}L^{H}Q^{\prime}}}} & (42) \\ {\overset{\_}{\beta} = \sqrt{\frac{E_{s}}{{{\hat{F}\left( {\text{:},1} \right)}}_{2}^{2} + {\sigma_{v}^{2}\;{{\hat{F}\left( {\text{:},{2\text{:}K}} \right)}}_{F}^{2}}}}} & (43) \end{matrix}$

This alternative manner of determining the feedback matrix and the pre-coding matrix determines the elements of these two matrices in a step by step manner, starting with the data stream that is to be pre-coded last. As a consequence this method provides more accurate results than the first method of calculating the precoding matrix F and feed-back matrix B. This improvement in accuracy, however, comes at the cost of increased computational complexity.

To compare the performance of the embodiments to other methods of determining receiver beam forming vectors, another method where each receiver's beam forming vector is the dominant left singular vector of its channel was evaluated. This method is referred to as ‘SVD’ in the figures. Also evaluated is an antenna selection method in which each receiver selects one channel, which has maximal norm among its channels. This method is referred to as ‘AS’ in the figures.

The performance of MMSE THP for the MISO case and ZF THP with receiver beam forming discussed in Min Huang et. al, (incorporated by reference in its entirety above) are also presented in the figures for comparison. The precoding and feedback matrices of the preferred embodiment are generated with user ordering based on equations (41) to (43).

The above methods were evaluated based on a simulation for a transmitter having four transmit antennae and communication with four receivers, each equipped with four receiver antennae. The simulations assume that the channel impulse response taps are independent and identically distributed (i.i.d.) with zero-mean Gaussian distribution. The simulation is further based on QPSK/16QAM modulation.

FIG. 5 illustrates the simulation results for QPSK modulation for receivers having four receiver antennae. It can be seen that the performances of the system using the receiver beam forming method of the preferred embodiment are considerably better than performance of the MMSE THP system for the MISO case (labelled 4×1 in FIG. 5) and that of the Zero-forcing THP proposed by Min Huang et. al. As can be seen from FIG. 3, the performance of the system operating based on the Cholesky based method of determining receiver beam forming vectors of the preferred embodiment provides considerable reductions in the bit error rate when compared to the systems operating based on the other methods.

FIG. 6 illustrates the simulation results for 16QAM modulation. It can again be seen that the system operating using the receiver beam forming vectors determined by the Cholesky based method of the preferred embodiment outperforms the comparative systems.

FIG. 7 illustrates simulation results for 16QAM modulation based on the assumption that the four users are selected from among eight potential users. For the ‘SVD’ and ‘AS’ schemes typical semi-orthogonal user selection is used. It will be appreciated from FIG. 5 that the performance of the proposed system using Cholesky-based user selection and receiver beam forming vector determination provides improvements over the ‘SVD’ and ‘AS’ schemes.

It will be appreciated that the above description of the embodiments is by way of example only. The person skilled in the art would readily appreciate modifications to this method. The method could, for example, be also applied to other MMSE-based precoding schemes based on Cholesky decomposition, such as vector perturbation precoding schemes, where the receiver beam forming vectors are generated firstly based on the Cholesky decomposition according to the method of a preferred embodiment and then vector perturbation precoding is conducted. One example for vector precoding can be found in “Minimum mean square error vector precoding”, Proc. 16^(th) IEEE Symp. on Personal, Indoor and Mobile Radio Communications, 2005 by D. Schmidt, M. Joham, and W. Utschick, the entirety of which is incorporated herein by this reference.

While certain embodiments have been described, the embodiments have been presented by way of example only, an area not intended to limit the scope of the inventions. Indeed, the novel methods and systems described herein may be embodied in a variety of other forms; furthermore, various omissions, substitutions and changes in the form of the methods and systems described herein may be made without departing from the spirit of the inventions. The accompanying claims and their equivalents are intended to cover such forms or modifications as would fall within the scope and spirit of the inventions. 

The invention claimed is:
 1. A method of determining receiver beam forming vectors for a plurality of receivers in a MIMO system, wherein each receiver comprises a plurality of receiver antennae and wherein the MIMO system comprises a plurality of transmitter antennae, the method comprising: quantifying the characteristics of the channels of the MIMO system in a channel property matrix Ψ including interference between individual channels, performing Cholesky decomposition of the channel property matrix row-by-row into a lower triangular matrix with unit elements along its diagonal, a diagonal matrix and the Hermitian transpose of the lower triangular matrix, comprising determining, when performing a step of the Cholesky decomposition for a row, a receiver beam forming vector associated with the row that maximises the component of the diagonal matrix in the row; wherein the matrix quantifying the characteristics of said channel is: ${H_{i}H_{i}^{H}} + {\gamma\; I_{N_{r}}} - {\sum\limits_{j = 1}^{i - 1}{s_{ij}s_{ij}^{H}q_{j}}}$ for i rows/steps of the Cholesky decomposition for i>0, wherein H_(i) is the channel matrix of the receiver associated with a current row/step of the Cholesky decomposition, H_(i) ^(H) is the Hermitian transpose of H_(i), γ is the expectation of noise that will be received I_(Nr) is an N_(r) ×N_(r) identity matrix, wherein N_(r) is the number of receiver antennae of the receiver, s_(ij) is $s_{ij} = \left\{ \begin{matrix} \frac{H_{i}{\overset{\Cap}{h}}_{i}^{H}}{q_{1}} & {{{if}\mspace{14mu} j} = 1} \\ \frac{\left( {{H_{i}{\overset{\Cap}{h}}_{j}^{H}} - {\sum\limits_{n = 1}^{j - 1}{s_{i\; n}l_{jn}^{*}q_{n}}}} \right)}{q_{j}} & {{{if}\mspace{14mu} j} > 1} \end{matrix} \right.$ ĥ_(j) ^(H) is the Hermitian transpose of r_(j)H_(j), wherein r_(j) is the receiver beam forming vector associated with the j-th element of the diagonal matrix, H_(j) is the channel matrix associated with the j-th element of the diagonal matrix, q_(i) is: $q_{i} = {\Psi_{ii} - {\sum\limits_{k}^{i - 1}{l_{ik}l_{ik}^{*}q_{k}}}}$ wherein: ${l_{ij} = {\frac{1}{q_{j}}\left( {\Psi_{ij} - {\sum\limits_{n}^{j - 1}{l_{i\; n}l_{jn}^{*}q_{n}}}} \right)}},{{{for}\mspace{14mu} i} > j}$ l_(ij) = 1, for  i = j l_(ij) = 0, for  i < j $\Psi_{ii} = {{{{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{i}^{H}} + {\gamma\mspace{14mu}{and}{\mspace{11mu}\;}\Psi_{ij}}} = {{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{j}^{H}}}$ ${\overset{\Cap}{h}}_{i} = {r_{i}H_{i}}$ s_(ij) ^(H) is the Hermitian transpose of s_(ij), l_(jn)* is the conjugate of l_(jn).
 2. A method according to claim 1, wherein maximising said component of the diagonal matrix comprises determining a left singular vector of the matrix quantifying the characteristics of the channel associated with an element of the diagonal matrix, wherein said left singular vector of the matrix is the left singular vector associated with the maximal singular value of the matrix.
 3. A method according to claim 1, wherein the matrix quantifying the characteristics of said channel quantifies channel gain, channel noise and said channel's interference with channels of receivers for which a beam forming vector had previously been determined.
 4. A method according to any of claim 1, wherein the matrices provided by the Cholesky decomposition of the channel property matrix Ψ are: $\begin{matrix} {\Psi = {LQL}^{H}} \\ {= {\begin{pmatrix} 1 & 0 & \ldots & 0 \\ l_{21} & 1 & \ldots & 0 \\ \ldots & \ldots & \ldots & \ldots \\ l_{K\; 1} & l_{K\; 2} & \ldots & 1 \end{pmatrix}\begin{pmatrix} q_{1} & 0 & \ldots & 0 \\ 0 & q_{3} & \ldots & 0 \\ \ldots & \ldots & \ldots & \ldots \\ 0 & 0 & \ldots & q_{K} \end{pmatrix}\begin{pmatrix} 1 & l_{21}^{*} & \ldots & l_{K\; 1}^{*} \\ 0 & 1 & \ldots & l_{K\; 2}^{*} \\ \ldots & \ldots & \ldots & \ldots \\ 0 & 0 & \ldots & 1 \end{pmatrix}}} \end{matrix}$ wherein L is the lower triangular matrix and has unit elements along its diagonal, Q is the diagonal matrix and L^(H) is the Hermitian transpose of L, wherein q₁ to q_(K) are the elements of the diagonal matrix defined as: $q_{i} = {\Psi_{ii} - {\sum\limits_{k}^{i - 1}{l_{ik}l_{ik}^{*}q_{k}}}}$ wherein: ${l_{ij} = {\frac{1}{q_{j}}\left( {\Psi_{ij} - {\sum\limits_{n}^{j - 1}{l_{i\; n}l_{jn}^{*}q_{n}}}} \right)}},{{{for}\mspace{14mu} i} > j}$ $\Psi_{ii} = {{{{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{i}^{H}} + {\gamma\mspace{14mu}{and}{\mspace{11mu}\;}\Psi_{ij}}} = {{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{j}^{H}}}$ ${\overset{\Cap}{h}}_{i} = {r_{i}H_{i}}$ wherein l_(ik)* is the conjugate of l_(ik), H_(i) is the channel matrix of the receiver associated with a current row/step of the Cholesky decomposition, γ is the expectation of noise that will be received and r_(i) is the receiver beam forming vector associated with the current row/step of the Cholesky decomposition.
 5. A method of determining a receiver beam forming vector for a receiver of a plurality of receivers in a MIMO system, wherein each receiver comprises plural receiver antennae, the method comprising: forming a matrix quantifying the characteristics of a channel of the receiver and interference said receiver experiences from other receivers for which a receiver beam forming vector had previously been determined, wherein the matrix quantifying the characteristics of said channel is: ${H_{i}H_{i}^{H}} + {\gamma\; I_{N_{r}}} - {\sum\limits_{j = 1}^{i - 1}{s_{ij}s_{ij}^{H}q_{j}}}$ for i rows/steps of the Cholesky decomposition for i>0, wherein H_(i) is the channel matrix of the receiver associated with a current row/step of the Cholesky decomposition, H_(i) ^(H) is the Hermitian transpose of H_(i), γ is the expectation of noise that will be received, I_(Nr) is an N_(r)×N_(r) identity matrix, wherein N_(r) is the number of receiver antennae of the receiver, s_(ij) is $s_{ij} = \left\{ \begin{matrix} \frac{H_{i}{\overset{\Cap}{h}}_{i}^{H}}{q_{1}} & {{{if}\mspace{14mu} j} = 1} \\ \frac{\left( {{H_{i}{\overset{\Cap}{h}}_{j}^{H}} - {\sum\limits_{n = 1}^{j - 1}{s_{i\; n}l_{jn}^{*}q_{n}}}} \right)}{q_{j}} & {{{if}\mspace{14mu} j} > 1} \end{matrix} \right.$ ĥ_(j) ^(H) is the Hermitian transpose of r_(j)H_(j), wherein r_(j) is the receiver beam forming vector associated with the j-th element of the diagonal matrix, H_(j) is the channel matrix associated with the j-th element of the diagonal matrix, q_(i) is: $q_{i} = {\Psi_{ii} - {\sum\limits_{k}^{i - 1}{l_{ik}l_{ik}^{*}q_{k}}}}$ wherein: ${l_{ij} = {\frac{1}{q_{j}}\left( {\Psi_{ij} - {\sum\limits_{n}^{j - 1}{l_{i\; n}l_{jn}^{*}q_{n}}}} \right)}},{{{for}\mspace{14mu} i} > j}$ l_(ij) = 1, for  i = j l_(ij) = 0, for  i < j $\Psi_{ii} = {{{{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{i}^{H}} + {\gamma\mspace{14mu}{and}{\mspace{11mu}\;}\Psi_{ij}}} = {{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{j}^{H}}}$ ${\overset{\Cap}{h}}_{i} = {r_{i}H_{i}}$ s_(ij) ^(H) is the Hermitian transpose of s_(ij), l_(jn)* is the conjugate of l_(jn); and determining the receiver beam forming vector including determining a left singular vector of the matrix that corresponds to the maximal singular value of the matrix.
 6. A method according to claim 1, further comprising determining a precoding and/or a feedback matrix for a transmitter using Tomlinson-Harashima precoding after two or more of the beam forming vectors have been determined.
 7. A device for determining receiver beam forming vectors for a plurality of receivers in a MIMO system, wherein each receiver comprises a plurality of receiver antennae and wherein the MIMO system comprises a plurality of transmitter antennae, the device comprising: a quantifier for quantifying in a channel property matrix the characteristics of receiver channels of the MIMO system, including interference between individual channels, a processor arranged to perform Cholesky decomposition of the channel property matrix in a row-by-row fashion into a lower triangular matrix with unit elements along its diagonal, a diagonal matrix and the Hermitian transpose of the lower triangular matrix, wherein the processor is further arranged to determine, when performing a step of the Cholesky decomposition for a row, a receiver beam forming vector associated with the row that maximises the component of the diagonal matrix in the row; wherein the matrix quantifying the characteristics of said channel is: ${H_{i}H_{i}^{H}} + {\gamma\; I_{N_{r}}} - {\sum\limits_{j = 1}^{i - 1}{s_{ij}s_{ij}^{H}q_{j}}}$ for i rows/steps of the Cholesky decomposition for i>0, wherein H_(i) is the channel matrix of the receiver associated with a current row/step of the Cholesky decomposition, H_(i) ^(H) is the Hermitian transpose of H_(i), γ is the expectation of noise that will be received, I_(Nr) is an N_(r)×N_(r) identity matrix, wherein N_(r) is the number of receiver antennae of the receiver, s_(ij) is $s_{ij} = \left\{ \begin{matrix} \frac{H_{i}{\overset{\Cap}{h}}_{j}^{H}}{q_{1}} & {{{if}\mspace{14mu} j} = 1} \\ \frac{\left( {{H_{i}{\overset{\Cap}{h}}_{j}^{H}} - {\sum\limits_{n = 1}^{j - 1}{s_{in}l_{jn}^{*}q_{n}}}} \right)}{q_{j}} & {{{if}\mspace{14mu} j} > 1} \end{matrix} \right.$ ĥ_(j) ^(H) is the Hermitian transpose of r_(j)H_(j), wherein r_(j) is the receiver beam forming vector associated with the j-th element of the diagonal matrix, H_(j) is the channel matrix associated with the j-th element of the diagonal matrix, q_(i) is: $q_{i} = {\Psi_{ii} - {\sum\limits_{k}^{i - 1}{l_{ik}l_{ik}^{*}q_{k}}}}$ wherein: ${l_{ij} = {\frac{1}{q_{j}}\left( {\Psi_{ij} - {\sum\limits_{n}^{j - 1}{l_{in}l_{jn}^{*}q_{n}}}} \right)}},\mspace{14mu}{{{for}\mspace{14mu} i} > j}$ l_(ij) = 1,  for  i = j l_(ij) = 0,  for  i < j $\Psi_{ii} = {{{{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{i}^{H}} + {\gamma\mspace{14mu}{and}\mspace{14mu}\Psi_{ij}}} = {{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{j}^{H}}}$ ĥ_(i) = r_(i)H_(i) s_(ij) ^(H) is the Hermitian transpose of s_(ij), l_(jn)* is the conjugate of l_(jn).
 8. A device according to claim 7, wherein the processor is arranged to maximise the component of the diagonal matrix by determining a left singular vector of the matrix quantifying the characteristics of the channel associated with an element of the diagonal matrix, wherein said left singular vector of the matrix is the left singular vector associated with the maximal singular value of the matrix.
 9. A device for determining a receiver beam forming vector for a receiver of a plurality of receivers in a MIMO system, wherein each receiver comprises plural receiver antennae, the device comprising: a quantifier for quantifying in a matrix the characteristics of a channel of the receiver and interference said receiver experiences from other receivers for which a receiver beam forming vector had previously been determined, wherein the matrix quantifying the characteristics of the said channel is: ${H_{i}H_{i}^{H}} + {\gamma\; I_{N_{r}}} - {\sum\limits_{j = 1}^{i - 1}{s_{ij}s_{ij}^{H}q_{j}}}$ for i rows/steps of the Cholesky decomposition for i>0, wherein H_(i) is the channel matrix of the receiver associated with a current row/step of the Cholesky decomposition, H_(i) ^(H) is the Hermitian transpose of H_(i), γ is the expectation of noise that will be received, I_(Nr) is an N_(r)×N_(r) identity matrix, wherein N_(r) is the number of receiver antennae of the receiver, s_(ij) is $s_{ij} = \left\{ \begin{matrix} \frac{H_{i}{\overset{\Cap}{h}}_{i}^{H}}{q_{1}} & {{{if}\mspace{14mu} j} = 1} \\ \frac{\left( {{H_{i}{\overset{\Cap}{h}}_{j}^{H}} - {\sum\limits_{n = 1}^{j - 1}{s_{in}l_{jn}^{*}q_{n}}}} \right)}{q_{j}} & {{{if}\mspace{14mu} j} > 1} \end{matrix} \right.$ ĥ_(j) ^(H) is the Hermitian transpose of r_(j)H_(j), wherein r_(j) is the receiver beam forming vector associated with the j-th element of the diagonal matrix, H_(j) is the channel matrix associated with the j-th element of the diagonal matrix, q_(i) is: $q_{i} = {\Psi_{ii} - {\sum\limits_{k}^{i - 1}{l_{ik}l_{ik}^{*}q_{k}}}}$ wherein: ${l_{ij} = {\frac{1}{q_{j}}\left( {\Psi_{ij} - {\sum\limits_{n}^{j - 1}{l_{in}l_{jn}^{*}q_{n}}}} \right)}},\mspace{14mu}{{{for}\mspace{14mu} i} > j}$ l_(ij) = 1,  for  i = j l_(ij) = 0,  for  i < j $\Psi_{ii} = {{{{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{i}^{H}} + {\gamma\mspace{14mu}{and}\mspace{14mu}\Psi_{ij}}} = {{\overset{\Cap}{h}}_{i}{\overset{\Cap}{h}}_{j}^{H}}}$ ${\overset{\Cap}{h}}_{i} = {r_{i}H_{i}}$ s_(ij) ^(H) is the Hermitian transpose of s_(ij), l_(jn)* is the conjugate of l_(jn); and a processor arranged to determine the receiver beam forming vector, wherein the processor is arranged to determine a left singular vector of the matrix that corresponds to the maximal singular value of the matrix. a processor arranged to perform Cholesky decomposition of the received matrix in a row-by-row fashion into a lower triangular matrix with unit elements along its diagonal, a diagonal matrix and the Hermitian transpose of the lower triangular matrix, wherein the processor is further arranged to determine, when performing a step of the Cholesky decomposition for a row, a receiver beam forming vector associated with the row that maximises the component of the diagonal matrix in the row.
 10. A device according to claim 7, wherein the processor is further arranged to determine a precoding matrix and/or a feedback matrix for a transmitter using Tomlinson-Harashima precoding after two or more of the receiver beam forming vectors have been determined.
 11. A non-transitory computer readable medium comprising computer executable instructions which, when executed by a computer, cause said computer to perform a method in accordance with claim
 1. 